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Abstract 

We propose a simple continuum model to interpret the shearing 
motion of dense, dry and cohesion- less granular media. Compressibil- 
ity, dilatancy and Coulomb- like friction are the three basic ingredients. 
The granular stress is split into a rate-dependent part representing 
the rebound-less impacts between grains and a rate-independent part 
associated with long-lived contacts. Because we consider stationary 
flows only, the grain compaction and the grain velocity are the two 
main variables. The predicted velocity and compaction profiles are in 
apparent agreement with the experimental or numerical results con- 
cerning free-surface shear flows as well as confined shear flows 

1 Introduction 

The mechanical behaviour of a flowing granular material depends strongly 
on the grain compaction. While dense granular media usually exhibit rela- 
tively slow motions with predominance of friction, less dense ones are usually 
found in vigorous motions with predominance of two-particles collisions. The 
collision-dominated regime is well described by kinetic theory, with the con- 
cepts of granular temperature and inelastic collisions. On the contrary, the 
current description of dense granular flows is not so fully satisfactory. It must 
be understood that we are not questioning the description by soil mechanics 
of quasi-static and highly stressed granular materials, but the description of 
flows with relatively low stress levels encountered, for example, in avalanches 



down an inclined plane. Several recent works (see e.g. [H 121 El) concluded 
rather pessimistically about the possibility of describing dense granular flows 
within the realm of continuum mechanics. In fact, the experimental observa- 
tion that most dense flows display a typical thickness of a few grain diameters 
must not be a factor of pessimism. We know from several examples in sus- 
pension mechanics that the continuum approach can cope with high velocity 
gradients in one direction, provided one has some statistical homogeneity in 
the other two directions. This situation is exactly the one met in sheared 
granular media, provided we discard transient effects and focus on the final 
stationary state. Once the continuum description is accepted, the number of 
relevant field variables must be decided. There is no doubt that the grain 
velocity is relevant and it is not less clear that the grain compaction is also a 
pertinent variable. In fact the widely used assumption of an incompressible 
medium is not tenable. It contradicts the dilatancy concept and, as will be 
seen below, the transport coefficients of a dense granular medium display 
enormous variations with only tiny modifications of the compaction. Our 
aim is thus to propose a model for dense and stationary shear flows in which 
the grain compaction and the grain velocity are the two fundamental vari- 
ables. One could also suggest the fiuctuational kinetic energy of the grains 
(the granular "temperature") as a third variable. However, since we limit 
our analysis to stationary shears, the granular temperature is no longer an 
independent variable. The role of the embedding fluid will be neglected ev- 
erywhere, and for these " dry" granular media, the main issue is to propose a 
constitutive relation for the granular stress. To compare with previous works 
on dense flows, we can say we adopt a phenomenological description some- 
what similar to that proposed two decades ago by Savage ^ and by Johnson 
and Jackson j3]. Like these authors, we introduce a stress tensor split into a 
frictional and a coUisional contribution. However, the coUisional contribution 
is concerned with rebound-less impacts characteristic of high grain concen- 
tration, and is free of any restitution coefficient. Our constitutive relation 
for the particulate stress has a form somewhat similar to that proposed by 
Ancey and Evesque jH] , the main differences concerning the role of the grain 
compaction and a more detailed expression of the granular pressure. Our 
model also shares some common features with the model proposed by Boc- 
quet et al. and by Louge 0, but instead of extending the kinetic theory 
approach to large compaction, we prefer here to develop a model specifically 
devoted to dense media. Moreover, the model we propose is quite simple in 
so far as it denies any special role to the compaction gradient P| and avoids 
the non-locality concept [lOj. 

Discarding two-particles collisions and any restitution coefficient means 
that our model is restricted to volume fractions in the range between (pm 
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and (pM- The maximum grain compaction (pM corresponds to the highest 
possible random packing (with (pM — 0.80 for two-dimensional flows and 
(j)M — 0.65 for three-dimensional ones) while (pm is the smallest compaction 
compatible with the existence of a continuous network of contacts between 
grains. As suggested by Azanza PP, one can define <pm as the minimum 
compaction for which the two-particle distribution function exhibits some 
swelling at a distance of two diameters. With this definition, 0^ — 0.70 for 
two-dimensional flows while (pm — 0.50 for three-dimensional ones. 

A phenomenological order parameter description of granular media was 
recently proposed pT|. We acknowledge this approach looks efficient in de- 
scribing a large number of phenomena observed in dense flows, but it suffers 
from two serious drawbacks. The first one is the huge number of possible 
candidates for the order parameter. For example, in the model to be pre- 
sented below, the reduced compaction (0 — 4>m)/{(pM — <Pm) could play this 
role. A second difficulty concerns the relevance of the Ginzburg-Landau re- 
laxation equation for the order parameter. Here we consider the standard 
conservation of mass and momentum only, without any extra equation. 

The description of stationary free-surface shear flows is given in section 
121 while that of conflned shear flows is presented in section El Section |31 
compares the model predictions with experimental and (or) numerical data. 
The flnal section insists on the limitations and necessary improvements of 
the proposed model, which must be considered as a minimal one. 



2 Free surface shear-flows 

As a prototype of shear flow with free surface, we consider the gravity-induced 
chute (over a heap or an inclined plate, see Figure with an angle 6 relative 
to the horizontal plane. The mean grain velocity is parallel to the x-axis, 
V = Ve^, while V and the solid fraction depend only on z, the distance 
to the free surface. The granular stress tensor is noted r and the equations 
of motion are: 

= + <Ppgsin{e) , = + ^pgcos{e) (1) 

where p is the constant mass per unit volume of the grain material and g is 
the acceleration of gravity. 

For dense granular media, the granular stress is a consequence of long- 
lived contacts and bounce-less impacts between grains. Long-lived contacts 
result from compressive forces acting towards the boundaries of the gran- 
ular medium. In the geometry considered, they take part in Tzz since z is 
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the direction of main compression. Whether gravity is responsible for com- 
pressive forces or not, we choose to scale the compressive stress with pgD 
where D is the grain size. The compressive stresses are related to the grain 
volume fraction as pgDF{(j)), where dF/dcj) is the non-dimensional compress- 
ibility of the granular medium. In free-surface shear flows, gravity is the only 
source of compaction and the magnitude of the compressive stress will also 
depend on 9. It is clear that the compressive role of gravity is maximum 
when the compression axis z is vertical while this role vanishes when gravity 
is orthogonal to it. Consequently, the general form of the gravity- induced 
compressive stress is pgDF{(l))f{6) with /(O) = 1 and f{n/2)= 0. The ex- 
act expression of f{6) is not important because, as will soon be seen, the 
stationary flows exist in a very limited range of 9 only. One of the simplest 
function of 9 which meets the above requirement is cos(6'), and we assume 
henceforth that the contribution of long-lived contacts to r^^ can be written 
in the form pgDF{(j))cos{9). To this gravity-induced contact stress must be 
added a rate-dependent impact stress. On purely dimensional grounds, this 
second contribution cannot be but Bagnold-like and the full normal stress 
finally appears in the form: 



where /XAr(0) represents the compaction-dependent intensity of the normal 
stress induced by the shear rate. Concerning the shear stress of the flowing 
granular medium, we assume it is made of a Coulomb-like contribution with 
a friction coefficient p{4>) completed by a Bagnold-like contribution involving 
a coefficient Pt{<P) representing the compaction-dependent intensity of the 
shear stress induced by the shear rate 



The model expressions 1^ and Q contain four functions of the grain com- 
paction. Before giving them some explicit (and tentative) expressions, let us 
comment on their expected general behaviour. These four functions are char- 
acteristic of the dense regime and have a meaning in the range (pm < < 4>m 
only. We expect F, pT and pn to become infinite when = 0m, because 
no motion nor extra compaction is expected above the maximum random 
packing. We also expect F and piy to vanish for = 0^, because the nor- 
mal stresses must vanish for the most tenuous contact network. Concerning 
the friction coefficient p, it is the only coefficient which remains finite when 
= 0M and it presumably increases for smaller compactions. In short. 






(3) 
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the three scalars F, fijsf and ht are strongly increasing functions of the com- 
paction, while /i has a much smoother behaviour. 

Since we neglect the role of the embedding fluid, the granular stress must 
vanish at the free surface and consequently t^z = tan(^)r^2 everywhere. In 
this case, when solving the equations of motion (P) with the model expres- 
sions (j2I) and 0, one arrives at a compaction profile and a velocity profile 
which are solution of: 



dz 9f F 1 W 



and 



9(/) L l-(/i]v/AtT)(tan(e)-/i) J 

DV^^dV f F(sin(^) - /icos(0)) '^^'^ 



(5) 

At the free-surface the solid fraction is ipm (remember we limit the description 
to the dense regime). According to (jl} the solid fraction increases towards 
its maximum value 0m over a depth which scales with the grain diameter but 
depends on 6 if UnIiJ'T is different from zero. Hence UnIiJ'T represents the 
relative magnitude of Reynold's dilatancy. Concerning the velocity profile, its 
characteristic value scales like {gDy/'^ and according to © its solution exists 
for any angle 6 verifying the inequality < tan(6') < -|-/iT(0)//^7v(0)- 
For certain values of 6 this inequality is possibly satisfied in a part only of 
the full range 4>m ^ 4> ^ 0a/- 

It is obviously not evident to deduce four functions of the compaction 
from the rather scarce experimental or numerical results on stationary shear 
flows. We assume henceforth that fi and f^r/f^'N are independent of the grain 
compaction. Then, a stationary solution is possible in a well-defined angle 
range 6'™„ < < 9max, with tan{9min) = and tan{9max) = /i + fir/fJ'N- 
To obtain more quantitative results we consider separately the chute over a 
heap from that over an inclined plane. 



2.1 Heap flows 

In the heap case, provided and i^n I I^^t are independant of the solid fraction, 
one can deduce from (j3)) and (jSl) the total granular flux flowing down the heap 
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the grain velocity Vheapi^) at the free surface 



and the relative velocity profile 



1/2 



fpM (f ^ dF i 



Vheapi^) r't'M f F \^^'^ dF ' 



C<f>M f 
■J (I'm \ 



Since the free-surface velocity and the flux are expected to have finite values 
for Omin < 9 < Omax, the two fuuctious F{(f)) and /ir(0) must be such as to 
guarantee the convergence of the above integrals. In this case Vheapi^) and 
Qheap are function of 9 with numerical prefactors depending on one's peculiar 
choice for F and /ir- In what follows we adopt the simple expressions 

F = FoLog — — and ht = /"to -7 t ■ (9) 

\ (pM - (P J \ 4>M -(P J 

The same expression for F was already proposed by Savage [U HH] and leads 
to an exponential-like volume fraction profile: 

with 

^ ^ 0M(l-^(tan(^^)-/x)) ^ ^ 

L{9) represents the typical thickness of the layer flowing down the heap, and 
it increases from L{9min) = to infinity when 9 = 9max- The relative 

velocity profile is also exponential-like for ~ 2 (see figure |21) but displays 

a Bagnold-like region of inverse concavity for ~ 0.2 (see figure E}. In 
fact, the numerical solution can be fitted by the analytical expression 



^/.eap(O) V 0A/-0m 7 \ 1 + (^ - l)e--AW 

With 0m = 0.5 and = 0.65 the total flux flowing down the heap is 

Qheap _ 14.^0^' (sin(^^)-/iCOs(e))V2 



(12) 



/iTo (i - B^(t^n{9) - fi) 



The dependence on of L and Qheap are represented in Fig. IH with /i = 0.36 
and fiN/fJ'T = 4.7. 
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2.2 Chute on rough plates 

The flows over inclined rough plates are more difficult to handle because the 
constitutive equations Q and (jHl) hold only in the bulk of the dense granular 
medium and are likely to be modified close to the rough plate. Since the 
role of the plate rugosity is difficult to assess quantitatively, we discard the 
description of the "basal layer" close to the plate El and assume a slip 
velocity Vg at some distance S above the rough plate. Then we apply ^ and 
(jHI) to a layer of thickness h, so that the free surface is located at a distance 
h + 6 above the rough incline. In the slab of thickness h, the solid fraction 
increases from 0^ at the free surface to the value 4>heap{h) at a distance 6 
from the rough plate where the velocity is Vg. The total flux through the 
core region is now given by: 

Qplate F{(j)h(,ap{h)) Vg 



l-^(tan(^)-/i) 



+ 



^ (sin(^) - /icos(^))i/2 f'f'^^-vW / p3\y^ Qp ^^^^ 



(tan(6') — /i) 



A rough plate is likely to slow down the core region more efficiently than a 
heap would do and we expect Vs < Vheapih). As a consequence, Qpiate{h, 6) as 
given in (|T!?|l is not expected to exceed Qheapifi) given in (jH]). When forcing 
a flux Qplate to flow dowu a rough plane inclined at angle 9, two different 
situations are encountered: when Qpiate is larger than Qheap{9), the granular 
medium will rearrange so as to flow down over a heap of angle 6 + a with 
Qplate = Qheap{0 + «)■ This gives a possible explanation for the "immature 
sliding flows" that were observed in some experiments jSl E] . Due to the very 
large increase of Qheap with 6 (see figure (jH)) and because the experimental 
flux is limited to some maximal value, immature sliding flows were observed 
for small angles close to O^nm only. Conversely, when Qpiate is smaller than 
Qheap{0) the whole layer of thickness h is in motion with a velocity everywhere 
larger than Vg. Moreover, when h/L{d) ~ 0.2, the Bagnold-like velocity 
proffie (which could hardly be observed in heap flows, see fig. (jH))) is now 
invading the whole core region. In fact, when expressions (jH)) are taken for 
granted and h/L{9) ~ 0.2, the total flux (|T3|l has the special form 

5/2 

D^gD D^gD ' 5^™ V m J \D J ' ^^^^ 

When the role of the velocity slip can be neglected, the second contribution 
gives the h^^"^ scaling which seems to be corroborated by experiments [TI] . 
as well as numerical simulations (IHj . 



Qplate _ (pmVsh 2 ( (sin(6') - /icos(6')) \ ^'"^ / h\ ' 
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3 Confined shear fiow 



In the two-dimensional shear flows we will consider the pressure load exerted 
on the boundaries of the granular medium is supposed to be applied along 
direction which is thus the direction of main compression. Because gravity 
plays a minor role concerning the compressive forces, the constitutive relation 
for is simply (compare with ^) 

— j +pgDF{cj>), (15) 

whatever the angle 6 between the z axis and gravity. The flow is along the 
X axis and the constitutive relation for the shear stress r^z is still given by 
((HI), without any change as compared to the free-surface case. 



3.1 Plane shear flow 

As a flrst type of conflned shear flow, we consider the planar shear of an 
inflnite horizontal granular layer bounded by two plates separated by a flxed 
distance h. The pressure load and the gravity are both oriented along the 
direction z and the flow is along direction x (see flgureE]). The equations of 
motion result in a constant shear stress S and a variable normal stress: 

= S and Tzz{z) = P(0) + pg f 0(0^^ 

where -P(O) is the pressure load exerted on the upper plate z = Q {z = h 
stands for the lower plate). We will distinguish the situation without and 
with gravity, the flrst case corresponding to numerical simulations and the 
second one to experiments. 



3.1.1 Without gravity 

In this case the normal stress is also a constant P all over the granular layer 
and the constitutive equations (fT^ and Q give: 

— j =P-pgDF{4>) (16) 

Depending on the sign of 5* — p{(f))P, we will have a static or a moving 
medium. In the static case the pressure load is noted Pq and the shear is 
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such that S < /x(0o)-Po where 0o is the constant compaction of the medium 
related to the pressure load through Pq = pgDF^ipQ). In the dynamic case, 
the compaction is still 0o (mass conservation) while the shear S is now larger 
than /i(0o)-Po- The velocity gradient is constant: 



dz J yUr(0o) + yu(0o)yUAf(0o) 

Due to dilatancy effects the pressure load exerted on the plates is necessarily 
larger than in the static case, following: 

P{S) = Po' ^-M0o)Po 



/i(0o) + /iT(0o)//WAf(0o) 

As a consequence the effective friction coefficient is a function (/)o and P: 

S /j^ ^ , /iT(0o) I, PdDj-,,,.. 
P (J-nWo) P 

3.1.2 With gravity 

In this case the normal stress increases in the downward direction so that 
the constitutive equation (jH)) results in 

pDVt(0) f ^ J =S- m(0)P(O) - fi{<P)pg 

It is then clear that the gravity-induced extra compaction possibly induces 
shear localization because the right-hand side can have a different sign in 
different parts of the flow. To simplify this issue we will now introduce the 
same assumptions we have previously used in the free-surface shear flows, 
namely that /i and pr/fJ'N do not depend on while and are given 

by (jni). We will flrst describe the static case before considering grain motions. 
Because the compaction on the upper plate is nesessarily different in the static 
and the dynamic cases, we deflne -Po(O) as the pressure load exerted on the 
upper plate when the granular medium is motionless and (j)o{z) as the static 
compaction proflle. As long as 5* < fiPo{0), the granular slab is motionless, 
the compaction 0o(O) at the upper plate satisfies Po(0) = pgDF{(f)Q{0)) and 
the compaction profile is: 

(j)o(z) = 7 ^^'^ . with Lq = —^D 



■>M 
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When the granular medium is flowing, the compaction profile (f){z) displays 
larger gradients and becomes 



with L 



_ 1^ p-./L 1 + 



^+U(0) ^ 



where (/)(0) is the new compaction at the upper plate. Since mass conservation 
requires 

j-h 

'(t){z) - (l)o{z)]dz = 0, 







it is clear that L < Lq results in 0(0) < 0o(O) and > (po{h). The 

compaction of the moving medium is thus reduced at the upper plate as 
compared to its static value while it is enhanced at the lower plate. 
The velocity profile is then deduced from the compaction profile 

1 + /i- 



J g \dz J ^t{4>) 

where S* is the dimensionless shear Let us define the volume fraction 6* 

pgD ^ 

such that S* = ^F{(j)*). It is clear that 0* > 0o(O) because S > /iPo(O). How- 
ever the above equation implies that motion exists for compactions less than 
0* only. This condition leads to check the self- consistency relation 0(z) < 0* 
for < z < h. This condition is automatically satisfied in the upper part 
of the flow since 0(0) < 0o(O) < 0*. But it may be not in the lower part, 
thus leading to a shear localization. This localization phenomenon is here de- 
pending on S* and h/L. Fig (jHI) shows the compaction and velocity profiles 
for two different values of 0*(S'*) investigating the two different situations 
depending whether localization occurs or not. 



3.2 Vertical chute flows 

A second type of confined shear flow is the chute between two vertical plates 
(see figure [7j). The compaction is due to a pressure load P exerted on the 
two plates along direction z. The flow and the gravity are oriented along 
direction x. The equations of motion result in a constant normal stress and 
a variable shear stress, by contrast to the preceding case: 

= P and r^^ = pg HOdL 
Jo 

where z = corresponds to the symmetry plane located between the two 
plates at which the shear stress vanishes. The constitutive relation (jH)) implies 
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pDV(0)f^) =pg md^-K<P)P- 

Either the right-hand side is everywhere negative (due to a very high pressure 
load) and the medium is motionless or there is a central region of the flow 
in which the shear stress does not exceed /iP and consequently where the 
strain rate vanishes. In this plug flow regime the solid fraction is a constant 
0* related to the pressure load as P = pgDF{(f)*). The thickness z* of the 
plug flow depends on 0* (hence on the pressure load) 

= ^P(r). 

Close to the vertical plates, there is a shear layer where the velocity decrease 
to Vui dependent on the plate roughness. In this parietal shear layer, the 
constitutive equations (fTHjl and (jHI) imply: 

D\ '/'dV^_ / P(0*)-P(0) \ 
gj dz { p^{4>) J ^ > 

and 



dz a 



To obtain more definite results we again consider the assumptions already 
made for gravity-driven and plane shear fiows, namely that fi and Pn/pt 
are independant of the solid fraction while P(0) and Pt{4') are given by 0. 
Then, the compaction profile in the shear layer z* < z < z^ is: 

^(z) = ^^^^ (19) 



where L* is the typical shear layer thickness: 

L* PtFq 



D pj^cj) 



M 



For the fiow to be dense up to the vertical plates, the wall compaction 
must be larger than 0m and the shear layer thickness is 

z -z* (^-^^ 
= Log ^ 
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As a consequence, the distance between the two plates is a function of 
0* (hence of P) and of (hence of the plate roughness). Concerning the 
velocity, it increases from a value at the wall to a value V^iug in the central 
part. The computed relative velocity field is represented in Fig (jH)) together 
with the fit 



4 Comparison with experimental and (or) 
numerical data 

4.1 Plane shear flow 

Neglecting the influence of gravity (as was done in most numerical simula- 
tions) our model leads to a uniform solid fraction and to a uniform velocity 
gradient, in conformity with results presented in figures 5b and 5c of [12] for 
the dense flow regime. When gravity is taken into account, a shear local- 
ization is possible, depending on the magnitude of the pressure load as well 
as on the thickness of the granular layer. Unfortunately, we are unaware of 
experimental or numerical data with which the predictions of Fig. © could 
be tested. 

4.2 Vertical chute flow 

The uniform solid fraction and the uniform velocity in the core region are 
correctly reproduced by the model. Concerning the sheared regions closed to 
the vertical boundaries, the relative velocity profile (^0]) and the compaction 
profile (fTTHl are quite similar to results presented in fig. 7b and 7c of [ 16^ . 

4.3 Heap flow 

The solid fraction profile ([111)) and the velocity profile (fT^ are quite close 
to those represented in fig. 9b and 9c of PH] and in fig. 9a and 9b of JZj. In 
particular, the velocity profile displays a Bagnold-like profile very close to 
the free-surface i^z < 0.2L{6)), a quasi-linear profile in the central part of 
the flow (0.2 < z/L{6) < 2) and finally an exponential tail for the deepest 
parts of the flow, observed in ^Hj. At variance with confined flows for which 
the shear was localized in boundary layers with thickness of the order of a 
few grain diameters, heap flows are characterized by a thickness L{6) of a 
few grain diameters when 6 is slightly larger than 6min but which increases 




(20) 
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to quite large values when 6 is close to 9max- A similar unlimited increase of 
the grain flux Qheap is observed for 9 close to 9max, as seen in fig. HI Such a 
behaviour is difficult to observe experimentally due to the limited values of 
Qheap that can be achieved in usual laboratory devices. According to (jB]) and 
(|TT|l . our model predicts L ~ {QheapY^^ when 9 is not close to 9min, a result 
slightly different from the scaling L ~ {QheapY^'^ suggested by fig.9j of [T^ . 

4.4 Rough inclined planes 

We explained the appearance of the so-called "immature sliding flows": they 
develop when the imposed flux Qpiate over a plate with inclination 9 is larger 
than the flux Qheap{9) which would fall down a heap with the same slope. 
Since Qpiate is experimentally limited to some maximum value, immature 
flows are observed for 9 close to 9min only. When 9 comes close to 9max, the 
thickness h of the granular layer flowing over the rough incline is smaller than 
the thickness of the grain layer which would flow down a heap with similar 
slope. And when h is less than 0.2L{9), the Bagnold-like velocity profile is 
invading the whole flowing layer, with the h^^'^ scaling for the flux Qpiate as a 
direct consequence (provided the first contribution to (fT^ is negligible). The 
main drawback of our model is its inability to explain the quantity hstop{9) 
introduced by Pouliquen |15j and which was confirmed in numerical simula- 
tions Pl]. The first reason is that we assumed the friction coefficient /i to be 
independent of the solid fraction. As a consequence 9min is a constant and 
hstop vanishes as soon as > 9min- A second reason is the possible inade- 
quacy of our model close to the rough incline. In this basal or frictional layer 
jHl IE!, the particle rotation plays an important role, the grain stress tensor 
is possibly non-symmetric and the solid fraction has a perturbed profile. All 
these phenomena would require a specific modelling. In fact the explanation 
of hstop{9) proposed by Mills et al. involves constitutive relations which 
are different close to the boundaries from those holding in the bulk. 

4.5 Annular shear 

This special kind of shear flow was not considered here because to describe 
it, we would need to give a constitutive equation for the Txx component of 
the granular stress, besides those for t^z and t^z- This will be done in future 
work. 
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5 Conclusion 



We proposed a model for dense granular flows which considers the solid 
fraction as the main microstructural parameter. The granular stress is par- 
titioned in a way similar to that proposed by Savage O E]. One of the 
distinctive features is a completely explicit expression for the contact stress 
which involves a function F{(j)) of the solid fraction. The solid fraction pro- 
file mainly depends on the compressibility dF/dip while the velocity profile 
is bound to F{(f))/ ^t{(P) where fiTi<P) is somehow analoguous to the effective 
viscosity used in the kinetic theory approach of dilute granular flows [HIH]. In 
principle the complete model contains two more functions of the solid fraction 
and /iiv(0)) but we strived to show that not so bad predictions could 
be obtained after assuming the friction coefficient /i and the dilatancy ratio 
A^tv/a^t to be independent of the solid fraction. Obviously, these are simplify- 
ing assumptions which can be released and improved. We also checked that 
the tentative (and simple) expressions for and iit{.<P) Isd to sound 
predictions. Needless to say that these expressions also can be improved. 
The main drawback of constitutive equations (j2)) and is their possible 
failure in a thin layer close to rough boundaries. Their main advantage is 
to contain all the ingredients necessary to interpret the shear-localization 
phenomenon, and to be able to explain the quite different velocity profiles 
appearing in the stationary shear flows of dense granular materials. 
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Figure 1: Gravity-induced shear flow with free surface (over a heap or an 
inclined plate with an angle 6 relative to the horizontal plane). 
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Figure 2: Reduced velocity profile Vheap/VheapiS^) versus the adimensional 
distance z/L{0) to the free surface. The dashed curve represents the approx- 
imate expression (fT^. The reduced compaction profile (0 — 0m) /(0m — 4>m) 
is plotted as well. 
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Figure 3: Zoom of figure |21 showing the Bagnold hke region for z/L{9) < 0.2. 
The dashed curve represents approximation ()12|) . 
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Figure 5: Planar confined shear flow of an inflnite horizontal granular layer 
bounded by two plates separated by a fixed distance h. The pressure load 
P and the gravity are both oriented along the direction z, the shear stress S 
and the flow are along direction x). 




Figure 6: Velocity proflle V{z) and compaction proflle <f){z), left large S, right 
small S. 
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Figure 7: Vertical chute flow between two rough plates. 
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K;)/(V^pi«5 - K,) versus the 
adimensional distance ^^^^ to the center. The arrow is in the increasing 

pgDFo' "^^^ dashed curve represents the approximate expression 
(f^ . The reduced compaction profiles 4>/(f>M with (pw/(pM = 0.7 are plotted 
as well. 



Figure 8: Reduced velocity profiles (y{z) 
adimensio 
values of 
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